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Abstract 

We examine the jp photoproduction and the hadronic 77 total cross sections by means of a QCD 
eikonal model with a dynamical infrared mass scale. In this model, where the dynamical gluon mass 
is the natural regulator for the tree level gluon-gluon scattering, the 'jp and 77 total cross sections 
are derived from the pp and pp forward scattering amplitudes assuming vector meson dominance 
and the additive quark model. We show that the validity of the cross section factorization relation 
(Jpp/a^p = a^p/a-y-y is fulfilled depending on the Monte Carlo model used to unfold the hadronic 
77 cross section data, and we discuss in detail the case of (7(77 hadrons) data with W^^ > 10 
GeV unfolded by the Monte Carlo generators PYTHIA and PHOJET. The data seems to favor a 
mild dependence with the energy of the probability {Phad) that the photon interacts as a hadron. 
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I. INTRODUCTION 



The hadronic nature of the photon is currently a subject of intense theoretical and exper- 
imental interest. The central question for hadronic interactions of real photons is whether 
they behave the same as hadrons or whether the additional hard contributions to the total 
cross sections of photon-induced interactions lead to a faster rise of the total 77 and 7p cross 
sections function of the energy. 

The increase of hadron-hadron total cross sections was theoreticallypredicted many years 
ago and this prediction has been accurately verified by experiment ^|. The measurements 
of the total photoproduction cross sections a"'^ at HERA N, \4i and the measurements of the 

nn 

total hadronic cross sections a'^'^ at LEP p, |6|| have also established the increase of photon- 
hadron and photon-photon total cross sections. Early modeling of these cross sections within 
Regge theory shows a energy dependence similar to the ones of nucleon-nucleon 0]. This 
universal behavior, appropriately scaled in order to take into account the differences between 
hadrons and the photon, can be understood as follows: at high center-of-mass energies the 
total photoproduction a'^^ and the total hadronic cross section a'^'^ for the production of 
hadrons in the interaction of one and two real photons, respectively, are expected to be 
dominated by interactions where the photon has fluctuated into a hadronic state. Therefore 
measuring the energy dependence of photon-induced processes should improve our under- 
standing of the hadronic nature of the photon as well as the universal high energy behavior 
of total hadronic cross sections. 

However the comparison of the experimental data and the theoretical prediction, as dis- 
cussed some time ago j^, may also present some subtleties depending on the Monte Carlo 
model used to analyze the data. Despite the similar energy dependence of hadron-hadron, 
photon-hadron and photon-photon total cross sections, there are significant differences be- 
tween photon induced cross section measurements and those of the purely hadronic processes. 
The latter comes from a counting rate AN{t), the number of counts/sec/ At in a narrow 
interval around the four-momentum transfer squared t, extrapolated to the value AN{t = 0) 
and normalized by an appropriate factor (for colliding beams this normalization factor is 
the luminosity L) . The 77 cross sections are extracted from a measurement of hadron pro- 
duction in e^e^ processes and are strongly dependent upon the acceptance corrections to 
be employed. These corrections are in turn sensitive to the Monte Carlo models used in 



2 



the simulation of the different components of an event, and this general procedure produces 
uncertainties in the determination of a"''^ ^. This clearly implies that any phenomenological 
analysis has to take properly into account the discrepancies among a"'"' data obtained from 
different Monte Carlo generators. 

A global description of photon-photon, photon-hadron and hadron-hadron total cross 
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. In these models 



sections is possible by means of QCD- inspired eikonal models j9, 
the increase of the total cross sections is associated with semihard scatterings of partons in 
the hadrons. The energy dependence of the cross sections is driven mainly by gluon-gluon 
scattering processes, since g{x, Q"^) ^ Q^) at small-x values. Nevertheless, the gluon- 
gluon subprocess cross section is potentially divergent at small transferred momenta, and the 
usual procedure to regulate this behavior is the introduction of a purely ad hoc parameter 



scale 



ating 

m 



separating^ the perturbative from the non-perturbative QCD region, like an infrared mass 
^ |lO, 
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15|. These arbitrary 



,114, 

or a cut-off at low transverse momentum px 
parameters are adjusted in order to obtain the best fits to the experimental data. 

Recently we introduced a QCD-inspired eikonal model where the ad hoc infrared mass 
scale was substituted by a dynamical gluon mass one Q|. One of the central advantages 
of the model is that it gives a precise physical meaning for the quoted infrared scale. Fur- 
thermore, since the behavior of the running coupling constant is constrained by the value of 



17|, the model also has a smaller number of parameters than similar 



dynamical gluon mass 
models. 

In this work we perform a detailed study of the 7p and 77 scattering in the framework of 



the QCD eikonal model with a gluon dynamical mass of Ref. jl6!]. We address the question 
of derivating 7p and 77 total cross sections from the pp and pp forward amplitudes, exploring 
the process underlying the extraction of 77 from the e^e~ data. The paper is organized 
as follows: in the next section we briefly review the QCD-inspired eikonal model with a 
dynamical gluon mass. Sec. Ill contains the details of 77 and 7p cross-sections calculations. 
The results are presented in the Sec. IV, where we provide fits for the cross sections, discuss 
the factorization theorem, app/a^p = cr^p/a^^, and investigate a possible mild dependence of 
the probability that the photon interacts as a hadron {Phad) increases with the energy. The 
results are compared to the data handled with the PYTHIA and PHOJET Monte Carlo 
generators. The conclusions are drawn in Sec. V. 
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II. A QCD-INSPIRED EIKONAL MODEL WITH A DYNAMICAL SCALE 



The high-energy hadron-hadron cross sections are calculated using a recently developed 
QCD-inspired eikonal model, where the onset of the dominance of gluons in the interaction 
of high-energy hadrons is managed by the dynamical gluon mass scale The model 

satisfies analyticity and unitarity constraints, and the hadron-hadron elastic scattering is 
due to the diffractive shadow of inelastic events. In the eikonal formalism the total cross 
section is given by 

atot{s) = 4n f bdb[l-e-^'^''''^cosx^{b,s)], (1) 

Jo 

where s is the square of the total center-of-mass energy, b is the impact parameter, and 
x{b, s) = Xr{^i s) + ixA^i s) is a complex eikonal function. Following the Ref. [l^, we write 
the even eikonal as the sum of gluon-gluon, quark-gluon, and quark-quark contributions: 

X'^b.s) = Xqq{b, S) + Xqgib, S) + Xggib, S) 

= i[(Jgq{s)W{b; figg) + aqg{s)W{b] {s)W{b;ftgg)], (2) 

where Xpp{b, s) = x~^{b, s) ±x^{b, s). Here W{b; /i) is the overlap function at impact param- 
eter space and crjj(s) are the elementary subprocess cross sections of colliding quarks and 
gluons = q, g). The overlap function, normalized so that J (PbW{b; /x) = 1, is associated 
with the Fourier transform of a dipole form factor, 

W{b-f,) = ^{f,b)^K,{^,b), (3) 

where K^i^x) is the modified Bessel function of second kind. The odd eikonal s), that 
accounts for the difference between pp and pp channels, is parametrized as 

X-(6,s) = C-S^e-/^W^(6;/i-), (4) 

where rrig is the dynamical gluon mass and the parameters and /^^ are constants to be 
fitted. The factor S is defined as 



(5) 



with the dynamical coupling constant set at its frozen infrared value. The origin of the 



dynamical gluon mass and the frozen coupling constant can be traced back to the ear 
of Cornwall and the formal expressions of these quantities can be seen in Ref. 



work 



Since that the main contribution to the asymptotic behavior of hadron-hadron total 
cross sections comes from gluon-gluon semihard coUisions, the eikonal functions Xqq{b,s) 
and Xqg{b,s), needed to describe the lower-energy forward data, are simply parametrized 
with terms dictated by the Regge phenomenology: 

TTl 

X,,ib,s)=tJ:A^W{b;^i,,), (6) 



Xgg{b, S) =iE 



A' + B'\n{A^ 

ml 



^{b\ y/figgfigg), (7) 



where A, A', B', fiqq and fitting parameters. 

In the model the gluon eikonal term Xgg{b,s) (see the expression is written as 
Xgg{b,s) = (T^^^{s)W{b;figg), where 

^r(^) = C'f drFgg{r)a^^\s). (8) 

J irrij^ I s 

Here Fgg^r) is the convoluted structure function for pair gg, (^gg^{,s) is the subprocess 
cross section and C is a fitting parameter. In the above expression it is introduced the 
energy threshold s > 4m^ for the final state gluons, assuming that these are screened gluons 
jigj ]. The structure function -Fgg(r) is given by 

Fggir) = [g ^ g]{r) = f ^ oi^) 9 {£l , (9) 



where g{x) is the gluon distribution function, adopted as 



g(r.) = N, <i-^, (10) 

where J = 1 + e and Ng = 2^0(6 — e)(5 — e)...(l — e). The correct analyticity properties of the 
model amplitudes is ensured by substituting s — *• se"""^/^ throughout Eqs. (0), © and (jHl). 
Other details of the model can be seen in Ref. and the fitting parameters are shown in 
Table m in the case of a dynamical gluon mass rUg = 400 MeV, which is the preferred value 
for pp and pp scattering. 



III. PHOTON-PROTON AND PHOTON-PHOTON REACTIONS 



The total cross section atot, the ratio of the real to the imaginary part of the forward 
scattering amplitude p, and the nuclear slope B for proton-proton and antiproton-proton 
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collisions have been found to be well described by the QCD-inspired eikonal model discussed 
in the last section, which we shall refer simply as the DGM model. Hence we fix the input 
parameters to the eikonal model calculations by using the data on pp and pp reactions in 
order to make predictions for (7"'"^ and a"''^ using the same values of the parameters. These 
values are shown in the Table HI 

The fact that at high energies a photon can fluctuate itself into a hadronic state open many 
theoretical possibilities. For example, in the model of Block et al. the data on 7p total 
photoproduction and the total hadronic 77 cross sections can be derived from the pp and pp 
forward amplitudes using vector meson dominance (VMD) and the additive quark model. 
We verify that this procedure, valid for QCD inspired models with an arbitrary infrared scale, 
can also be applied to our eikonal model with a dynamical gluon mass being necessary only 
minimal changes displayed in the sequence. Thus, considering that a strongly interacting 
photon behaves as a system of two quarks (the additive quark model), the even amplitude 
for 7p scattering can be obtained after the substitutions aij | o"jj and /ij^ i^i 
the eikonal of Eq. (j21), resulting in the following expression: 



2 

X+ (6, s) = 



ggyS)W\ b; Xll^igg 



.11) 



2 

Note that in the work of Block et al. the above procedure is used to obtain directly 
the physical eikonal x^^i^js), what means that the odd part of the eikonal that describes 
pp and pp scattering is neglected, and x'^^{b,s) = x^p{b,s). This is a good approximation 
at high energies, where the total cross sections pp and pp are asymptotically expected to be 
the same. However, below -^i ~ 100 GeV the pp and pp total cross sections are different, 
as happens in other channels like n~p and n^p as well as K~p and K^p. Therefore it is 
reasonable to expect the same behavior for the channels jp and jp. To account for these 
possible differences at low energy in the case of a'^^ and o"'^^ we apply the prescription 
dictated by the additive quark model also in the odd eikonal of Eq. (jH), 

X^,{b, s) = lc-^^^ e-/^ W L W|/i- I , (12) 



and the full physical eikonal, combining the even and odd parts, will be finally given by 
Xfp = X^p{b,s)±x~p{b,s). 

Assuming vector meson dominance (VMD), the eikonalized 'jp total cross section can be 
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written as 



bdbpp[l - e-^?'(^'") cosx7(&, ^)]' 



(13) 



where P^J^^ is the probabihty that the photon interacts as a hadron. In the simplest VMD 
formulation this probability is expected to be of 0{aem)- 



had 



P, 



had 



V=p,ijj,i 



249 



(14) 



where p, uo and are vector mesons. However, there are expected contributions to Phad 
other than p, u;, (f). For example, in the eikonalized minijet model analysis of a"'"' of the Ref. 
jllj l a good fit to the photoproduction data is obtained with the value Phad = 1/204, which 
includes a non-VMD contribution of (9(20%). In this case there are strong contributions 
of heavier vector mesons and continuum states. Moreover, the probability Phad may also 
depend on the energy, which is a possibility that we will explore in this paper. 
To extend the model to the 77 channel we just perform the substitutions o"jj 



|(Tjj and 



\lJiij in the even part of the eikonal (Eq. Q): 
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(15) 



hence the additive quark model and the VMD are again assumed to be valid, and the 
calculation leads to the following eikonalized total 77 hadronic cross section 



a 



77/ 
tot\ 



^had 



bdb [1 - e-^r('''^) cosx^^(&, s)] 



(16) 



where 



p77 
had 



had' 



(17) 



With the eikonal parameters of the DGM model fixed by the pp and pp data (Table H]), 
we have performed all calculations of photoproduction and photon-photon scattering. Apart 
from the equations that we have displayed above, we will introduce two possible modifica- 
tions when we present our results in the next section. One is related to a feasible energy 
dependence of Phad) and the other is concerned with the requirement of a normalization 
factor in the a"''^ calculation. 
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IV. RESULTS 



In our analysis the center-of-mass energy thresholds for the 7p and 77 channels are the 
same that one used in the study of the pp and pp channels: -y/i™" = W^™" = 10 GeV. We 
adopt a fitting procedure, assuming an interval ~ xLin corresponding, in the case of 
normal errors, to the projection of the x^ hypersurface containing 90% of probability. In the 
case of 1, 2 and 3 free parameters, this corresponds to the interval x"^ ~Xmin = 2.70, 4.61 and 



6.25, respective 
Group (PDG) 



/. We use the data set for a'^^ compiled and analyzed by the Particle Data 
I, with the statistic and systematic errors added in quadrature. Data on 
77 total hadronic cross sections are also available on the PDG compilation, but the 77 data 
quoted by the PDG are the average of the results obtained by unfolding the data with the 



PYTHIA j20[ and the PHOJET 2l| Monte Carlo generators. These event generators, based 
on a model by Schuler and Sjostrand |,22j (PYTHIA) and on the Dual Parton model j^l 
(PHOJET), are used to simulate photon-photon interactions and, in general, data unfolded 
with PYTHIA are higher than if unfolded with PHOJET. As a consequence we obtain totally 
different scenarios for the rise of the total hadronic 77 cross section at high energies. For 
this reason, we have deconvoluted the data on a'^'^ of the L3 and OPAL p collaborations, 
according to whether PYTHIA or PHOJET was used. Therefore, instead of using the PDG 
data for a^'^, we performed global fits considering separately the data obtained through the 
PYTHIA and PHOJET codes, defining two data sets as 

SET I: a^P and afyj, data (v^^^, W^^ > 10 GeV), 
SET II: a^P and aJ^HO data (v^^^, W^^ > 10 GeV), 

where (t'Pyt {^Vho) correspond to the data of 77 total hadronic cross section obtained via 
the PYTHIA (PHOJET) generator. The unfolded results are depicted in the Figure 1. This 
strategy of analysis is similar to the one adopted by Block and Kang, where a global test of 
factorization for nucleon-nucleon, 7p and 77 scattering was performed through real analytic 
amplitudes [24]. Their results have evidenced the critical role played by the Monte Carlo 
simulations used by L3 and OPAL group. 

Our first analysis (henceforth referred to as FIT I) consisted in the determination of the 
Phad values from fits to the data sets defined above. In the fit to the SET I we obtained 
p-^i = 220.85 ± 0.59; the x^lDOF for this fit was 3.29. The fit to the SET II has resulted 
in P-^^ = 223.70 ± 0.50, with x^/DOF = 1.94. These values are shown in Table II. The 
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and a'^'^ cross sections corresponding to these values of P^ad shown by the sohd curves in 
Figure 2. The results of the fit to the SET I are displayed in Figs. 2(a) and 2(b), whereas 
the results of the fit to the SET II are displayed in Figs. 2(c) and 2(d). In the FIT I the 
/DOF values are relatively high. However, as we already commented, the extraction of 
77 data from e^e~ processes is strongly sensitive to the Monte Carlo simulation, and this 
procedure could require a systematic overall-normalization factor in the experimental data. 
In fact, the cross section (7(77 hadrons) is determined from the differential cross section 
of the reaction e+e^ e"'"e~7*7* e'^e" hadrons by means of both the luminosity function 
fo. tKe pWo„ fl.x and the Kadro.c section .^IQ^ BBHy Q. 
However, the uncertainty in the cross section extrapolation to real photons {Qf = Q2 = ^) 
is estimated to be 5-7.5%, and this uncertainty is not included in the systematic error of 



the measurement [27[. In this case we can introduce a normalization factor in the a'^'^ data 
as performed in Ref. ^| or, what we have chosen to do, just multiply Eq. (fTH|l by a 
normalization factor to be fitted by the experiment: 



77 

'-'tot 



» = ATrNP^a, I hdh[l- e-^r(M cosxr(6, s)]. (18) 

J n 



Within this approach (henceforth referred to as FIT II) we carried out new fits to the 
SETS I and II described above. We obtained for the SET I (SET II) the values P^^^ = 
222.30 ± 1.18 (P;-^^ = 220.30 ± 1.71) and N = 1.086 ± 0.028 (iV = 0.936 ± 0.025), with 
X^/DOF = 2.05 {x^/DOF = 1.02). These values are shown in Table II. The jp and 77 
total cross sections corresponding to these values of and are shown by the dashed 
curves in Figure 2. The results of the fit to the SET I are displayed in Figs. 2(a) and 
2(b), whereas the results of the fit to the SET II are displayed in Figs. 2(c) and 2(d). We 
notice that the approach adopted in the FIT II leads to an improvement of the fits, since 
the x^/DOF values obtained through the FIT II are smaller than ones obtained through 
the FIT I. These results indicate the convenience of introducing a normalization factor in 
order to account for possible effects not included in the systematic error of a"'^ data. 

In the sequence we have also explored the possibility that Phad varies with the energy. 
Such possibility has been raised in the past 128], although there are arguments against it 
I29I 1. We have assumed a phenomenological expression for Phad, implying that it increases 
logarithmically with the square of the center of mass energy: Phad = a + b\n{s). Our 
intention was just to see if the quality of the fit improves in such situation, referred to 
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as FIT III. The results of the fit to the SET I were: a = (3.882 ± 0.122) x 10-^ b = 
(1.195 ± 0.229) X 10'^ and = 0.963 ± 0.037, with x^/DOF = 0.86. The fit results to the 
SET II were: a = (4. 175 ± 0.247) x IQ-^, b = (6.272 ± 4.775) x 10"^ and = 0.875 ± 0.052, 
with x^/DOF = 0.76. All these resuhs are shown in Table III and the total cross section 
curves are depicted in Figure 3, where Figs. 3(a) and 3(b) are related to the SET I whereas 
the Figs. 3(c) and 3(d) are related to the SET II. Considering the smaller /DOF values 
obtained when Phad increases with the energy (FIT III), we believe that this possibility 
cannot be completely neglected. 

In all the cases we see that the data of (tJ^yt above ^/s ~ 100 GeV can ha rdly be described 
by our model, but this is also a problem for other models in the literature [20| based on the 
factorization assumption 



Assuming the correctness of these models we may interpret the steeper rise of the a"''^ 
appointed by the PYTHIA as a possible signal of violation of the factorization relation. Oth- 
erwise this behavior indicates a possible failure of the PYTHIA generator at high-energies. 
The latter possibility is qualitatively supported by the results depicted in the Figures 3(c) 
and 3(d), that show that the shape and normalization of the PHOJET cross sections are in 
good agreement with the factorization relation (jl9p . 



V. CONCLUSIONS 



Our first aim and main result was to show that an eikonal QCD inspired model based on 
a nonperturbative infrared dynamical gluon mass scale with the help of vector meson 
dominance and the additive quark model, can successfully describe the data of the total 
photoproduction 7p and total hadronic 77 cross sections. The calculation is similar to the 
one of Block et al. [9], with the advantage that in the DGM model one has a physical 
meaning for the infrared scale, and the coupling constant is not a free parameter, turning 
out to be related to the dynamical gluon mass scale. 

We promoted some minor changes in the calculation, assuming that Phad-, the probability 
that a photon interacts as a hadron, has a logarithmic increase with s. This leads to a 
improvement in the fit, i. e. it lowers the x^/DOF values. There are pros [2^ and cons 
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29[ about this possibility. We just point out that the logarithmic increase of Phad with s is 
quite favored by the data. 

We verified that the data of cipyj. above y/s ~ 100 GeV can hardly be described by our 
model, but this is also a problem for other models in the literature. Assuming the correctness 
of these models we could say that the PHOJET generator is more appropriate to obtain the 
0"'*"*' data above ^/s ~ 100 GeV. Connected with this comment we can recall that recently 
Block and Kaidalov demonstrated that the factorization relation (fT^ does not depend on 
the assumption of an additive quark model, but more on the opacity of the eikonal being 
independent of the nature of the reaction ^jj. Hence, according to this result, we argue that 
it will be difficult for QCD-inspired models incorporating the total cross section factorization 
hypothesis to accommodate the (7(77 — >• hadrons) data with y^^^ > 100 GeV unfolded with 
the PYTHIA Monte Carlo generator. 
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TABLE I: Values of the parameters of the DGM model resulting from the global fit to the forward 
pp and pp data . The dynamical gluon mass scale was set to nig = 400 MeV. 



C (12.097±0.962)xl0-3 

Hgg [GeV] 0.7242±0.0172 

A 6.72±0.92 

Hgg [GeV] 1.0745±0.0405 

A' (4.491±0.179)xl0"3 

B' 1.08±0.14 

C- 3.17±0.35 

fi- [GeV] 0.6092±0.0884 





TABLE II: Fits adopting Phad 


constant. 






FIT I 

SET I SET II 


SET I 


FIT II 

SET II 


p-l 

had 

N 

xVdof 


220.85±0.59 223.70±0.50 
3.29 1.94 


222.30±1.18 
1.086±0.028 
2.05 


222.30±1.71 
0.936±0.025 
1.02 


FIGURES 


2(a), 2(b) 2(c), 2(d) 


2(a), 2(b) 


2(c), 2(d) 




TABLE III: Fits adopting Phad = 


= a + 61n(s). 






SET I 




SET II 


a 
b 

N 

xVdof 


(3.882±0.122)xl0"3 
(1.195±0.229)xl0-^ 
0.963±0.037 
0.86 




(4.175±0.247)xl0"3 
(6.272±4.775)xl0^5 
0.875±0.052 
0.76 


FIGURES 


3(a), 3(b) 




3(c), 3(d) 
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_ m L3, unfolded with PYTHIA 
- ■ OPAL, unfolded with PYTHIA 
1000 - o L3, unfolded with PHOJET 
' □ OPAL, unfolded with PHOJET 

800 - 



600 - 




200 - 

I I I I I I I I ^ I I 

10 10^ 

W„ (GeV) 

FIG. 1: cr'>'T data from the L3 and OPAL collaborations unfolded with the PYTHIA and PHOJET 
generators. 
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^0.24 
6 0.22 



(a) 




b^ooo : a, L3, unfolded with PYTHIA 
: n OPAL, unfolded with PYTHIA 



(b) 




Vs (GeV) 



Vs (GeV) 



(c) 




: o L3, unfolded with PHOJET 
: n OPAL, unfolded with PHOJET 



Vs (GeV) 




Vs (GeV) 



FIG. 2: a'^P and a'^'^ cross sections corresponding to the constant values of Pf^^^^ shown in Table II. 
The solid (dashed) curves correspond to the FIT I (II), whereas the results of the fit to the SET I 
are displayed in Figs. 2(a) and 2(b), and the results of the fit to the SET II are displayed in Figs. 
2(c) and 2(d). 
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^0.24 
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(a) 




b^ooo : a, L3, unfolded with PYTHIA 
: n OPAL, unfolded with PYTHIA 



(b) 




Vs (GeV) 



Vs (GeV) 



(c) 




: o L3, unfolded with PHOJET 
: n OPAL, unfolded with PHOJET 



Vs (GeV) 




Vs (GeV) 



FIG. 3: a'^'P and a'^'^ cross sections corresponding to the the case where varies with the energy 
as described in Table III. The curves of Figs. 3(a) and 3(b) are related to the fit labeled as SET I 
in Table III, whereas the Figs. 3(c) and 3(d) are related to the values of the SET II displayed in 
the same table. 
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